Individualized detection of TMPRSS2-ERG fusion status in prostate cancer: a rank-based qualitative transcriptome signature.

BACKGROUND: TMPRSS2-ERG (T2E) fusion is highly related to aggressive clinical features in prostate cancer (PC), which guides individual therapy. However, current fusion prediction tools lacked enough accuracy and biomarkers were unable to be applied to individuals across different platforms due to their quantitative nature. This study aims to identify a transcriptome signature to detect the T2E fusion status of PC at the individual level. METHODS: Based on 272 high-throughput mRNA expression profiles from the Sboner dataset, we developed a rank-based algorithm to identify a qualitative signature to detect T2E fusion in PC. The signature was validated in 1223 samples from three external datasets (Setlur, Clarissa, and TCGA). RESULTS: A signature, composed of five mRNAs coupled to ERG (five ERG-mRNA pairs, 5-ERG-mRPs), was developed to distinguish T2E fusion status in PC. 5-ERG-mRPs reached 84.56% accuracy in Sboner dataset, which was verified in Setlur dataset (n = 455, accuracy = 82.20%) and Clarissa dataset (n = 118, accuracy = 81.36%). Besides, for 495 samples from TCGA, two subtypes classified by 5-ERG-mRPs showed a higher level of significance in various T2E fusion features than subtypes obtained through current fusion prediction tools, such as STAR-Fusion. CONCLUSIONS: Overall, 5-ERG-mRPs can robustly detect T2E fusion in PC at the individual level, which can be used on any gene measurement platform without specific normalization procedures. Hence, 5-ERG-mRPs may serve as an auxiliary tool for PC patient management.

Up-regulated SPP1 increases the risk from IPF to lung cancer via activating the pro-tumor macrophages.

The incidence of lung cancer (LC) in Idiopathic Pulmonary Fibrosis (IPF) patients is more than twice that in non-IPF. This study aims to investigate IPF-to-LC pathogenesis and to develop a predictor for detecting IPF predisposing patients to LC. We conducted unsupervised clustering to detect high-risk subtypes from IPF to LC. Subsequently, we performed single-cell RNA-seq analysis to characterize high-risk IPF by examining the immune microenvironment. We identified 42 common immune function-related pathogenic genes between IPF and LC. We developed an LC risk classifier for IPF patients, comprising five genes: SPP1, MMP9, MMP12, FABP4, and IL1B. The five-gene classifier can successfully distinguish the high-risk population from IPF patients. High-risk IPF patients exhibited an immunosuppressive microenvironment with higher oncogene expression than low-risk patients. Single-cell analysis revealed that SPP1+ macrophages at the terminal of macrophages' developmental trajectory may promote the progression from IPF to LC. The strong crosstalk between SPP1+ macrophages and inflammation-related cancer-associated fibroblasts promoted the tumorigenic process in IPF. In vitro, assays showed that co-culturing macrophages overexpressing SPP1 with MRC-5 cells induced the transition of fibroblasts into cancer-associated fibroblasts. SPP1 produced by macrophages promoted epithelial-mesenchymal transition in alveolar epithelial cells via stimulating the upregulation of N-cadherin and Vimentin in MLE-12 cells. This study provided a novel method to identify the LC risk population from IPF, revealing the cellular interactions involved in the transition from IPF to LC. Our findings highlighted SPP1 as a critical driver in IPF progression, offering a potential target for therapy in fibrosis.

Multi-omics integration analysis unveils heterogeneity in breast cancer at the individual level.

Identifying robust breast cancer subtypes will help to reveal the cancer heterogeneity. However, previous breast cancer subtypes were based on population-level quantitative gene expression, which is affected by batch effects and cannot be applied to individuals. We detected differential gene expression, genomic, and epigenomic alterations to identify driver differential expression at the individual level. The individual driver differential expression reflected the breast cancer patients' heterogeneity and revealed four subtypes. Mesenchymal subtype as the most aggressive subtype harbored deletion and downregulated expression of genes in chromosome 11q23 region. Specifically, silencing of the SDHD gene in 11q23 promoted the invasion and migration of breast cancer cells in vitro by the epithelial-mesenchymal transition. The immunologically hot subtype displayed an immune-hot microenvironment, including high T-cell infiltration and upregulated PD-1 and CTLA4. Luminal and genomic-unstable subtypes showed opposite macrophage polarization, which may be regulated by the ligand-receptor pairs of CD99. The integration of multi-omics data at the individual level provides a powerful framework for elucidating the heterogeneity of breast cancer.

Synthetic viability induces resistance to immune checkpoint inhibitors in cancer cells.

BACKGROUND: Immune checkpoint inhibitors (ICI) have revolutionized the treatment for multiple cancers. However, most of patients encounter resistance. Synthetic viability (SV) between genes could induce resistance. In this study, we established SV signature to predict the efficacy of ICI treatment for melanoma. METHODS: We collected features and predicted SV gene pairs by random forest classifier. This work prioritized SV gene pairs based on CRISPR/Cas9 screens. SV gene pairs signature were constructed to predict the response to ICI for melanoma patients. RESULTS: This study predicted robust SV gene pairs based on 14 features. Filtered by CRISPR/Cas9 screens, we identified 1,861 SV gene pairs, which were also related with prognosis across multiple cancer types. Next, we constructed the six SV pairs signature to predict resistance to ICI for melanoma patients. This study applied the six SV pairs signature to divide melanoma patients into high-risk and low-risk. High-risk melanoma patients were associated with worse response after ICI treatment. Immune landscape analysis revealed that high-risk melanoma patients had lower natural killer cells and CD8+ T cells infiltration. CONCLUSIONS: In summary, the 14 features classifier accurately predicted robust SV gene pairs for cancer. The six SV pairs signature could predict resistance to ICI.

A genetic map of the chromatin regulators to drug response in cancer cells.

BACKGROUND: Diverse drug vulnerabilities owing to the Chromatin regulators (CRs) genetic interaction across various cancers, but the identification of CRs genetic interaction remains challenging. METHODS: In order to provide a global view of the CRs genetic interaction in cancer cells, we developed a method to identify potential drug response-related CRs genetic interactions for specific cancer types by integrating the screen of CRISPR-Cas9 and pharmacogenomic response datasets. RESULTS: Totally, 625 drug response-related CRs synthetic lethality (CSL) interactions and 288 CRs synthetic viability (CSV) interactions were detected. Systematically network analysis presented CRs genetic interactions have biological function relationship. Furthermore, we validated CRs genetic interactions induce multiple omics deregulation in The Cancer Genome Atlas. We revealed the colon adenocarcinoma patients (COAD) with mutations of a CRs set (EP300, MSH6, NSD2 and TRRAP) mediate a better survival with low expression of MAP2 and could benefit from taxnes. While the COAD patients carrying at least one of the CSV interactions in Vorinostat CSV module confer a poor prognosis and may be resistant to Vorinostat treatment. CONCLUSIONS: The CRs genetic interaction map provides a rich resource to investigate cancer-associated CRs genetic interaction and proposes a powerful strategy of biomarker discovery to guide the rational use of agents in cancer therapy.

LiqBioer: a manually curated database of cancer biomarkers in body fluid.

Cancer biomarkers are measurable indicators that play vital roles in clinical applications. Biomarkers in body fluids have gained considerable attention since the development of liquid biopsy, and their data volume is rapidly increasing. Nevertheless, current research lacks the compilation of published cancer body fluid biomarkers into a centralized and sustainable repository for researchers and clinicians, despite a handful of small-scale and specific data resources. To fulfill this purpose, we developed liquid biomarker (LiqBioer) containing 6231 manually curated records from 3447 studies, covering 3056 biomarkers and 74 types of cancer in 22 tissues. LiqBioer allows users to browse and download comprehensive information on body liquid biomarkers, including cancer types, source studies and clinical usage. As a comprehensive resource for body fluid biomarkers of cancer, LiqBioer is a powerful tool for researchers and clinicians to query and retrieve biomarkers in liquid biopsy. DATABASE URL: http://www.medsysbio.org:8080/LiqBioer.

Systematic exploration of the underlying mechanism of gemcitabine resistance in pancreatic adenocarcinoma.

Resistance to gemcitabine is the main challenge of chemotherapy for pancreatic ductal adenocarcinoma (PDAC). Hence, the development of a response signature to gemcitabine is essential for precision therapy of PDAC. However, existing quantitative signatures of gemcitabine are susceptible to batch effects and variations in sequencing platforms. Therefore, based on within-sample relative expression ordering of pairwise genes, we developed a transcriptome-based gemcitabine signature consisting of 28 gene pairs (28-GPS) that could predict response to gemcitabine for PDAC at the individual level. The 28-GPS was superior to previous quantitative signatures in terms of classification accuracy and prognostic performance. Resistant samples classified by 28-GPS showed poorer overall survival, higher genomic instability, lower immune infiltration, higher metabolic level and higher-fidelity DNA damage repair compared with sensitive samples. In addition, we found that gemcitabine combined with phosphoinositide 3-kinase (PI3K) inhibitor may be an alternative treatment strategy for PDAC. Single-cell analysis revealed that cancer cells in the same PDAC sample showed both the characteristics of sensitivity and resistance to gemcitabine, and the activation of the TGFß signalling pathway could promote progression of PDAC. In brief, 28-GPS could robustly determine whether PDAC is resistant or sensitive to gemcitabine, and may be an auxiliary tool for clinical treatment.

Upregulation of CXCL1 and LY9 contributes to BRCAness in ovarian cancer and mediates response to PARPi and immune checkpoint blockade.

BACKGROUND: Mutations in BRCA1 or BRCA2 (BRCA1/2) cause homologous recombination deficiency (HRD). Ovarian cancer (OvCa) patients harbouring HRD beyond BRCA1/2 mutation result in a state referred to as "BRCAness". OvCa with BRCAness could benefit from PARP inhibitors. This study aims to identify a signature to detect the BRCAness population at the transcriptome level. METHODS: We used a rank-based algorithm to develop a qualitative BRCAness signature for OvCa. Upregulation of CXCL1 with downregulation of SV2A and upregulation of LY9 with downregulation of CHRNB3 were constructed as the BRCAness signature (2 gene pairs, 2-GPS) for OvCa. RESULTS: OvCa samples that were classified as BRCAness by 2-GPS showed improved overall survival, progression-free survival and exhibited increased multi-omics alterations in homologous recombination genes and enhanced sensitivity to immune checkpoint blockade. BRCAness cells were sensitive to PARP inhibitors. By biological experiments, we validated SKOV3 cells and patients with HRD exhibited higher expression of CXCL1 than SV2A and higher expression of LY9 than CHRNB3 at mRNA level. Both SKOV3 and A2780 with HRD were sensitive to mitomycin C, cisplatin and olaparib. CONCLUSIONS: In conclusion, 2-GPS could robustly predict BRCAness OvCa at the individual level and extend the population who may benefit from PARP inhibitors.

Pan-cancer multi-omics analyses reveal crosstalk between the Hippo and immune signaling pathways in the tumor microenvironment.

The Hippo signaling pathway is critical for carcinogenesis. However, the roles of the Hippo signaling pathway in the tumor immune microenvironment have been rarely investigated. This study systematically analyzed the relationship between the Hippo signaling pathway and immune cell infiltration across 32 cancer types. Both bioinformatics analyses and biological experiments revealed that the downstream effector of Hippo signaling YAP1 might inhibit CD8+ T cell infiltration by upregulating the expression of the transcription factor CREB1 in uterine corpus endometrial carcinoma. In addition, esophageal carcinoma (ESCA) patients were classified into three subtypes based on the Hippo-immune gene panel. The subtypes of ESCA had distinct characteristics in immune cell infiltration, immune pathways, and prognosis. Thus, this study also reveals a new classification of the immune subtypes with prognostic characteristics in ESCA.

A transcriptional signature detects homologous recombination deficiency in pancreatic cancer at the individual level.

Pancreatic cancer (PC) with homologous recombination deficiency (HRD) has been reported to benefit from poly ADP-ribose polymerase (PARP) inhibitors. However, accurate identification of HRD status for PC patients from the transcriptional level is still a great challenge. Here, based on a relative expression ordering (REO)-based algorithm, we developed an HRD signature including 24 gene pairs (24-GPS) using PC transcriptional profiles from The Cancer Genome Atlas (TCGA). HRD samples classified by 24-GPS showed worse overall survival (p = 4.4E-3 for TCGA; p = 1.2E-3 for International Cancer Genome Consortium-Australia cohort; p = 6.4E-2 for GSE17891; p = 7.5E-2 for GSE57495) and higher HRD scores than non-HRD samples (p = 1.4E-4). HRD samples showed highly unstable genomic characteristics and also displayed HRD-related alterations at the epigenomic and proteomic levels. Moreover, HRD cell lines identified by 24-GPS tended to be sensitive to PARP inhibitors (p = 6.6E-2 for olaparib; p = 2.6E-3 for niraparib). Compared with the non-HRD group, the HRD group presented lower immune scores and CD4/CD8 T cell infiltration proportion. Interestingly, PC tumor cells with co-inhibition of PARP-related genes and ATR showed reduced survival ability. In conclusion, 24-GPS can robustly identify PC patients with HRD status at the individualized level.

Discovering a qualitative transcriptional signature of homologous recombination defectiveness for prostate cancer.

The discovery of homologous recombination deficiency (HRD) biomarkers in prostate cancer is important for patients who will benefit from poly ADP-ribose polymerase inhibitor (PARPi). Here, we developed a transcriptional homologous recombination defectiveness (HRDness) signature, comprising 16 gene pairs (16-GPS), for prostate cancer by a relative expression ordering (REO)-based discovery procedure. Subsequently, two newly subtypes classified by 16-GPS showed a higher significance level in various clinicopathological and HRD features than subtypes obtained by other methods, such as HRDetect. HRDness subtype also displayed more aggressive features and higher genomics scores than non-HRDness in three independent datasets. HRDness prostate cancer cells were more sensitive to PARPi than non-HRDness. Moreover, the HRDness samples showed distinct multi-omics characteristics related to homologous recombination repair function loss. Overall, the newly proposed qualitative signature can robustly determine the HRD status for prostate cancer at the personalized level, and especially be an auxiliary tool for PARPi treatment strategy.

Revealing biomarkers associated with PARP inhibitors based on genetic interactions in cancer genome.

Poly (ADPribose) polymerase inhibitors (PARPis) are clinically approved drugs designed according to the concept of synthetic lethality (SL) interaction. It is crucial to expand the scale of patients who can benefit from PARPis, and overcome drug resistance associated with it. Genetic interactions (GIs) include SL and synthetic viability (SV) that participate in drug response in cancer cells. Based on the hypothesis that mutated genes with SL or SV interactions with PARP1/2/3 are potential sensitive or resistant PARPis biomarkers, respectively, we developed a novel computational method to identify them. We analyzed fitness variation of cell lines to identify PARP1/2/3-related GIs according to CRISPR/Cas9 and RNA interference functional screens. Potential resistant/sensitive mutated genes were identified using pharmacogenomic datasets. We identified 41 candidate resistant and 130 candidate sensitive PARPi-response related genes, and observed that EGFR with gain-of-function mutation induced PARPi resistance, and predicted a combination therapy with PARP inhibitor (veliparib) and EGFR inhibitor (erlotinib) for lung cancer. We also revealed that a resistant gene set (TNN, PLEC, and TRIP12) in lower grade glioma and a sensitive gene set (BRCA2, TOP3A, and ASCC3) in ovarian cancer, which were associated with prognosis. Thus, cancer genome-derived GIs provide new insights for identifying PARPi biomarkers and a new avenue for precision therapeutics.

Individualized lncRNA differential expression profile reveals heterogeneity of breast cancer.

Long non-coding RNAs (lncRNAs) play key regulatory roles in breast cancer. However, population-level differential expression analysis methods disregard the heterogeneous expression of lncRNAs in individual patients. Therefore, we individualized lncRNA expression profiles for breast invasive carcinoma (BRCA) using the method of LncRNA Individualization (LncRIndiv). After evaluating the robustness of LncRIndiv, we constructed an individualized differentially expressed lncRNA (IDElncRNA) profile for BRCA and investigated the subtype-specific IDElncRNAs. The breast cancer subtype-specific IDElncRNA showed frequent co-occurrence with alterations of protein-coding genes, including mutations, copy number variation and differential methylation. We performed hierarchical clustering to subdivide TNBC and revealed mesenchymal subtype and immune subtype for TNBC. The TNBC immune subtype showed a better prognosis than the TNBC mesenchymal subtype. LncRNA PTOV1-AS1 was the top differentially expressed lncRNA in the mesenchymal subtype. And biological experiments validated that the upregulation of PTOV1-AS1 could downregulate TJP1 (ZO-1) and E-Cadherin, and upregulate Vimentin, which suggests PTOV1-AS1 may promote epithelial-mesenchymal transition and lead to migration and invasion of TNBC cells. The mesenchymal subtype showed a higher fraction of M2 macrophages, whereas the immune subtype was more associated with CD4 + T cells. The immune subtype is characterized by genomic instability and upregulation of immune checkpoint genes, thereby suggesting a potential response to immunosuppressive drugs. Last, drug response analysis revealed lncRNA ENSG00000230082 (PRRT3-AS1) is a potential resistance biomarker for paclitaxel in BRCA treatment. Our analysis highlights that IDElncRNAs can characterize inter-tumor heterogeneity in BRCA and the new TNBC subtypes indicate novel insights into TNBC immunotherapy.

Multilevel prioritization of gene regulators associated with consensus molecular subtypes of colorectal cancer.

Consensus molecular subtypes (CMSs) are emerging as critical factor for prognosis and treatment of colorectal cancer. Gene regulators, including chromatin regulator, RNA-binding protein and transcriptional factor, are critical modulators of cancer hallmark, yet little is known regarding the underlying functional mechanism in CMSs. Herein, we identified a core set of 235 functional gene regulators (FGRs) by integrating genome, epigenome, transcriptome and interactome of CMSs. FGRs exhibited significant multi-omics alterations and impacts on cell lines growth, as well as significantly enriched cancer driver genes and pathways. Moreover, common FGRs played different roles in the context of CMSs. In accordance with the immune characteristics of CMSs, we found that the anti-tumor immune pathways were mainly activated by FGRs (e.g. STAT1 and CREBBP) in CMS1, while inhibited by FGRs in CMS2-4. FGRs mediated aberrant expression of ligands, which bind to receptor on immune cells, and modulated tumor immune microenvironment of subtypes. Intriguingly, systematic exploration of datasets using genomic and transcriptome co-similarity reveals the coordinated manner in FGRs act in CMSs to orchestrate their pathways and patients' prognosis. Expression signatures of the FGRs revealed an optimized CMS classifier, which demonstrated 88% concordance with the gold-standard classifier, but avoiding the influence of sample composition. Overall, our integrative analysis identified FGRs to regulate core tumorigenic processes/pathways across CMSs.

Reference genome and annotation updates lead to contradictory prognostic predictions in gene expression signatures: a case study of resected stage I lung adenocarcinoma.

RNA-sequencing enables accurate and low-cost transcriptome-wide detection. However, expression estimates vary as reference genomes and gene annotations are updated, confounding existing expression-based prognostic signatures. Herein, prognostic 9-gene pair signature (GPS) was applied to 197 patients with stage I lung adenocarcinoma derived from previous and latest data from The Cancer Genome Atlas (TCGA) processed with different reference genomes and annotations. For 9-GPS, 6.6% of patients exhibited discordant risk classifications between the two TCGA versions. Similar results were observed for other prognostic signatures, including IRGPI, 15-gene and ORACLE. We found that conflicting annotations for gene length and overlap were the major cause of their discordant risk classification. Therefore, we constructed a prognostic 40-GPS based on stable genes across GENCODE v20-v30 and validated it using public data of 471 stage I samples (log-rank P < 0.0010). Risk classification was still stable in RNA-sequencing data processed with the newest GENCODE v32 versus GENCODE v20-v30. Specifically, 40-GPS could predict survival for 30 stage I samples with formalin-fixed paraffin-embedded tissues (log-rank P = 0.0177). In conclusion, this method overcomes the vulnerability of existing prognostic signatures due to reference genome and annotation updates. 40-GPS may offer individualized clinical applications due to its prognostic accuracy and classification stability.

An individualized transcriptional signature to predict the epithelial-mesenchymal transition based on relative expression ordering.

The epithelial-mesenchymal transition (EMT) process is involved in cancer cell metastasis and immune system activation. Hence, identification of gene expression signatures capable of predicting the EMT status of cancer cells is essential for development of therapeutic strategies. However, quantitative identification of EMT markers is limited by batch effects, the platform used, or normalization methods. We hypothesized that a set of EMT-related relative expression orderings are highly stable in epithelial samples yet are reversed in mesenchymal samples. To test this hypothesis, we analyzed transcriptome data for ovarian cancer cohorts from publicly available databases, to develop a qualitative 16-gene pair signature (16-GPS) that effectively distinguishes the mesenchymal from epithelial phenotype. Our method was superior to previous quantitative methods in terms of classification accuracy and applicability to individualized patients without requiring data normalization. Patients with mesenchymal-like ovarian cancer showed poorer overall survival compared to patients with epithelial-like ovarian cancer. Additionally, EMT score was positively correlated with expression of immune checkpoint genes and metastasis. We, therefore, established a robust EMT 16-GPS that is independent of detection platform, batch effects and individual variations, and which represents a qualitative signature for investigating the EMT and providing insights into immunotherapy for ovarian cancer patients.

Biased Influences of Low Tumor Purity on Mutation Detection in Cancer.

The non-cancerous components in tumor tissues, e.g., infiltrating stromal cells and immune cells, dilute tumor purity and might confound genomic mutation profile analyses and the identification of pathological biomarkers. It is necessary to systematically evaluate the influence of tumor purity. Here, using public gastric cancer samples from The Cancer Genome Atlas (TCGA), we firstly showed that numbers of mutation, separately called by four algorithms, were significant positively correlated with tumor purities (all p < 0.05, Spearman rank correlation). Similar results were also observed in other nine cancers from TCGA. Notably, the result was further confirmed by six in-house samples from two gastric cancer patients and five in-house samples from two colorectal cancer patients with different tumor purities. Furthermore, the metastasis mechanism of gastric cancer may be incorrectly characterized as numbers of mutation and tumor purities of 248 lymph node metastatic (N + M0) samples were both significantly lower than those of 121 non-metastatic (N0M0) samples (p < 0.05, Wilcoxon rank-sum test). Similar phenomena were also observed that tumor purities could confound the analysis of histological subtypes of cancer and the identification of microsatellite instability status (MSI) in both gastric and colon cancer. Finally, we suggested that the higher tumor purity, such as above 70%, rather than 60%, could be better to meet the requirement of mutation calling. In conclusion, the influence of tumor purity on the genomic mutation profile and pathological analyses should be fully considered in the further study.

RHPCG: a database of the Regulation of the Hippo Pathway in Cancer Genome.

The Hippo signaling pathway is a highly conserved pathway controlling organ size, cell proliferation, apoptosis and other biological functions. Recent studies have shown that Hippo signaling pathway also plays important roles in cancer initiation and progression. However, a database offering multi-omics analyses and visualization of Hippo pathway genes in cancer, as well as comprehensive Hippo regulatory relationships is still lacking. To fill this gap, we constructed the Regulation of the Hippo Pathway in Cancer Genome (RHPCG) database. Currently, RHPCG focuses on analyzing the 21 core Hippo-protein-encoding genes in over 10 000 patients across 33 TCGA (The Cancer Genome Atlas) cancer types at the levels of genomic, epigenomic and transcriptomic landscape. Concurrently, RHPCG provides in its motif section 11 regulatory motif types associated with 21 core Hippo pathway genes containing 180 miRNAs, 6182 lncRNAs, 728 circRNAs and 335 protein coding genes. Thus, RHPCG is a powerful tool that could help researchers understand gene alterations and regulatory mechanisms in the Hippo signaling pathway in cancer.

Genetic Interaction-Based Biomarkers Identification for Drug Resistance and Sensitivity in Cancer Cells.

Cancer cells generally harbor hundreds of alterations in the cancer genomes and act as crucial factors in the development and progression of cancer. Gene alterations in the cancer genome form genetic interactions, which affect the response of patients to drugs. We developed an algorithm that mines copy number alteration and whole-exome mutation profiles from The Cancer Genome Atlas (TCGA), as well as functional screen data generated to identify potential genetic interactions for specific cancer types. As a result, 4,529 synthetic viability (SV) interactions and 10,637 synthetic lethality (SL) interactions were detected. The pharmacogenomic datasets revealed that SV interactions induced drug resistance in cancer cells and that SL interactions mediated drug sensitivity in cancer cells. Deletions of HDAC1 and DVL1, both of which participate in the Notch signaling pathway, had an SV effect in cancer cells, and deletion of DVL1 induced resistance to HDAC1 inhibitors in cancer cells. In addition, patients with low expression of both HDAC1 and DVL1 had poor prognosis. Finally, by integrating current reported genetic interactions from other studies, the Cancer Genetic Interaction database (CGIdb) (http://www.medsysbio.org/CGIdb) was constructed, providing a convenient retrieval for genetic interactions in cancer.

Link synthetic lethality to drug sensitivity of cancer cells.

Synthetic lethal (SL) interactions occur when alterations in two genes lead to cell death but alteration in only one of them is not lethal. SL interactions provide a new strategy for molecular-targeted cancer therapy. Currently, there are few drugs targeting SL interactions that entered into clinical trials. Therefore, it is necessary to investigate the link between SL interactions and drug sensitivity of cancer cells systematically for drug development purpose. We identified SL interactions by integrating the high-throughput data from The Cancer Genome Atlas, small hairpin RNA data and genetic interactions of yeast. By integrating SL interactions from other studies, we tested whether the SL pairs that consist of drug target genes and the genes with genomic alterations are related with drug sensitivity of cancer cells. We found that only 6.26%∼34.61% of SL interactions showed the expected significant drug sensitivity using the pooled cancer cell line data from different tissues, but the proportion increased significantly to approximately 90% using the cancer cell line data for each specific tissue. From an independent pharmacogenomics data of 41 breast cancer cell lines, we found three SL interactions (ABL1-IFI16, ABL1-SLC50A1 and ABL1-SYT11) showed significantly better prognosis for the patients with both genes being altered than the patients with only one gene being altered, which partially supports the SL effect between the gene pairs. Our study not only provides a new way for unraveling the complex mechanisms of drug sensitivity but also suggests numerous potentially important drug targets for cancer therapy.